Article 


Polarized thermal emission from dustina 
galaxy at redshift 2.6 


https://doi.org/10.1038/s41586-023-06346-4 


Received: 17 February 2023 


Accepted: 20 June 2023 


Published online: 06 September 2023 


Open access 


® Check for updates 


J. E. Geach'™, E. Lopez-Rodriguez’, M. J. Doherty’, Jianhang Chen®, R. J. lvison®*>, 
G. J. Bendo’, S. Dye® & K. E. K. Coppin' 


Magnetic fields are fundamental to the evolution of galaxies, playing a key role inthe 
astrophysics of the interstellar medium and star formation. Large-scale ordered 
magnetic fields have been mapped in the Milky Way and nearby galaxies’”, but it 

is not known how early in the Universe such structures formed’. Here we report the 
detection of linearly polarized thermal emission from dust grains ina strongly lensed, 
intrinsically luminous galaxy that is forming stars at arate more than 1,000 times 
that of the Milky Way at redshift 2.6, within 2.5 Gyr of the Big Bang*>. The polarized 
emission arises from the alignment of dust grains with the local magnetic field®”. The 
median polarization fraction is of the order of 1%, similar to nearby spiral galaxies’. 
Our observations support the presence of a 5-kiloparsec-scale ordered magnetic field 
with a strength of around 500 uG or lower, oriented parallel to the molecular gas disk. 
This confirms that such structures can be rapidly formed in galaxies, early in cosmic 


history. 


We observed the lensed galaxy 9i09 (ref. 4) with the Atacama Large 
Millimeter/Submillimeter Array (ALMA) at a representative frequency 
of 242 GHz (equivalent to a wavelength of roughly 350 um in the rest- 
frame of the galaxy) to record the dust continuum emission averaged 
over a total bandwidth of 7.5 GHz. The set of XX, YY, XYand YX linear 
polarization parameters recorded in full polarization mode allow meas- 
urement of the Stokes parameters Q and U, yielding the total linearly 
polarized intensity, PI=./Q?+U? and position angle of polarized 
emission y= 0.5 arctan(U/Q). The root mean squared sensitivity of the 
observations is 0, = 47 and og S oy = 9 py beam”. In Fig. 1 we present 
image plane maps of the total intensity /, Stokes Qand U, and polarized 
intensity (PI). The polarization angle x is rotated by 90° to show the 
plane-of-the-sky magnetic (B) field orientation (xs). We measure an 
image plane integrated flux density of /= 62 mJy, integrated polariza- 
tion fraction of P= 0.6 + 0.1%, where P = PI// and B field orientation of 
Xp = (-0.7 +1.4)°. The mean of the distribution of polarization fractions 
and B field orientations is (P) = 0.6 + 0.3% and (xp) = (0.8 + 18.3)°, respec- 
tively. Note that the uncertainties are the dispersion of the distribution 
of individual measurements within the galaxy, not the accuracy inthe 
polarization measurement (Methods). 

Using the lens model derived from previous high-resolution milli- 
metre continuum emission and optical Hubble Space Telescope imag- 
ing, the source plane CO(4-3) emission, tracing the cold molecular gas 
reservoir, has been shown to be well modelled by a rotating disk of 
maximum radius 2.6 kpc, inclined by roughly 50° to the line of sight, 
witha position angle onthe sky of roughly 5° east of north*». With this 
model as a constraint, we explore what source plane B field configura- 
tions are consistent with the image plane polarization observations. 
The most likely source plane configuration is a large-scale ordered B 


+5° 


field oriented y, = 574, east of north with an extent matching that of 
the CO emission (Fig. 2). This result implies the presence of a5-kpc-scale 
galactic ordered magnetic field oriented parallel to the molecular 
gas-rich disk. Angular variations of y, across the galaxy present in the 
image plane maps, corresponding to scales of 600 pc in the source 
plane, can be explained by the low signal to noise ratio and beam effects 
(Methods). This result indicates that the introduction of arandom 
B field component with an angular variation of +5° in addition to the 
large-scale ordered B field is also consistent with the observations. 
Apresent, we lack the sensitivity and resolution to map the configura- 
tion of the B field at scales of roughly 100 pcin which structure related 
toturbulence can start to be resolved. The observed B field configura- 
tion parallel to the disk is consistent with the galactic B fields measured 
in local spiral galaxies observed at far-infrared and radio wavelengths’”. 
Note that our far-infrared polarimetric observations trace a density- 
weighted average B field in the cold and dense interstellar medium 
(ISM), rather than a volume average B field in the warm and diffuse ISM 
by radio polarimetric observations. 

The mean and integrated polarization fractions of 9i09 are consistent 
withthe P = 0.8% level measured in nearby spiral and starburst galaxies 
at wavelengths of 53-214 um (ref. 2). The observations presented here 
are sensitive to polarized emission beyond this range, pushing into the 
Rayleigh—Jeans tail of the thermal emission spectrum at a rest-frame 
wavelength of A,..= 350 pm. In recent models of diffuse interstellar 
dust, the polarization fraction, P, isindependent of wavelength across 
200-2,000 um, consistent with observations of Galactic dust emission. 
Observations of local starburst galaxies show that Ponly varies by 0.4% 
over the 50-150 um range, with an increase of up to roughly 1% towards 
214 um (ref. 2). We therefore conclude that 9i09 has a polarization 


'Centre for Astrophysics Research, School of Physics, Engineering and Computer Science, University of Hertfordshire, Hatfield, UK. "Kavli Institute for Particle Astrophysics and Cosmology, 
Stanford University, Stanford, CA, USA. “European Southern Observatory, Garching, Germany. “Department of Physics and Astronomy, Macquarie University, Sydney, New South Wales, 
Australia. School of Cosmic Physics, Dublin Institute for Advanced Studies, Dublin, Ireland. “Institute for Astronomy, Royal Observatory, University of Edinburgh, Edinburgh, UK. 7UK ALMA 
Regional Centre Node, Jodrell Bank Centre for Astrophysics, Department of Physics and Astronomy, The University of Manchester, Manchester, UK. °School of Physics and Astronomy, 


University of Nottingham, Nottingham, UK. “e-mail: j.geach@herts.ac.uk 


Nature | www.nature.com | 1 


Article 


a Stokes / b Stokes Q c Stokes U d 
0° 16’ 02” 0° 16’ 02” 0° 16’ 02” a 0° 16’ 02” 
` > 

Za. yee ao ao jg 
te 
Š - 
8 15’ 58” 15’ 58” 15 58” M J 15’ 58” 

eee 56” 56” eee 

3 4 

e f g h 
0° 16’ 02” 0° 16’ 02” 0° 16’ 02” E 0° 16’ 02” 
Ss 00” 00” oe 00” -$ 00” 
=] 
N 
3 
g 15’ 58” 15’ 58” 15’ 58” . 15’ 58” 
j . 

56” 56” 56” 56” 

O O O e- O O O O 
o »? 4 ©? ow? o? 9? N 
Nn » y » u Y » t D i D bi » K DA 
KY & & & 
gŭ 2 RA (J2000) 2 RA (J2000) 2 RA (J2000) 
< < < 
o g g LY 
0 2 4 6 8 10 -0.04 -0.02 0 0.02 0.04 -0.04 -0.02 0 0.02 0.04 0 0.02 0.04 0.06 0.08 0.10 


| (mJy per beam) Q (mJy per beam) 


Fig. 1 | The magnetic field orientation of the gravitationally lensed galaxy 
9i09 atz = 2.553. a-d, ALMA 242 GHz polarimetric observations of the Stokes 
I, Qand U parameters, andthe polarized intensity (PI). The synthetic beam of 
the observations (1.2” x 0.9”, 8 = 68°) is shown as the red ellipse, lower left. The 
B field orientation is indicated by white lines shown at the Nyquist sampling, 


level similar to local star-forming disks and starburst galaxies, witha 
key difference being the order-of-magnitude difference in gas mass 
and star-formation rate, with the disk of 9i09 being close to molecular 
gas dominated, contrasted with the f,,; = 10% gas fractions of local 
star-forming disks’. 

The large-scale ordered magnetic fields that exist in massive disk 
galaxies in the local Universe is thought to arise through the amplifica- 
tion of seed fields, and this has been predicted to occur on relatively 
short cosmological timescales, of order 1 Gyr (refs. 10-12). Weak seed 
fields (as low as B = 10°” G) could be formed in protogalaxies either 
through trapping of acosmological field, possibly primordial in nature, 
or through the battery effect following the onset of star formation” *. 
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Fig. 2 | Source plane configuration of the magnetic field and lensing model. 
a, Source plane intensity and field orientation. b, Lensed source plane image. 
c, Synthetic observations with the synthetic beam size (1.2” x 0.9”, 8 = 68°) 
indicated by the red ellipse. The B field orientation is indicated by white lines 
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with line lengths proportional to the polarization fraction. e-h, Synthetic 
polarimetric observations using a constant B field configuration in the source 
plane. Contours indicate signal to noise: for Stokes /, the contours increase as 
0, x 2345, For Stokes Qand Uand for PI, the contours start at 30 and increase in 
steps of lo. Dec., declination; RA, right ascension. 


Although turbulent gas motions in disks can reduce net polarization if 
they impart a strong turbulent component to the B field”, recent theo- 
retical models of the formation of galactic-scale magnetic fields invoke 
turbulence in the ISM as the origin of a ‘small-scale’ dynamo that can 
rapidly amplify the weak seed fields to uG levels”"*”. This small-scale 
dynamo is mainly driven by supernova explosions with coherence 
lengths of order 50-100 pc, but turbulence can be injected into the ISM 
on multiple scales through disk instabilities and feedback effects, includ- 
ing stellar winds and outflows driven by radiation pressure, supernova 
explosions and large-scale outflows from an active galactic nucleus. 
The average turbulent velocity component of the disk of 9i09, deter- 
mined from kinematic modelling of the CO emission, is o, = 70 kms * 


Synthetic observation 


(Pip) = 1.140.1% 
(pip) = 4.3 + 6.3° 


Offset (^) 


with lengths proportional to the polarization fraction. The median and root 
mean squared values of the polarization fraction and B field orientation are 
indicated in the images. The caustics in the source plane and image plane are 
shown as green and yellow lines, respectively. 


and the star-formation rate density exceeds 100 M, yr™ kpc” (ref. 5). 
The high dense gas fraction of the molecular reservoir—as traced by 
the ratio of CO(4-3)/C 1(1-O) emission—is also consistent with the 
injection of supersonic turbulence, which plays a key role in shaping 
the lognormal probability distribution function of the molecular gas 
density”. There is also tentative evidence of stellar feedback inaction 
through the broad lines of dense gas tracers’. Finally, one expects a 
high cosmic-ray flux density in the ISM of 9i09, commensurate with the 
high star-formation rate density, and this too could serve to amplify 
magnetic fields. Therefore, 9i09 probably has the conditions required 
to rapidly amplify any weak seed fields by means of the small-scale 
dynamo effect, with amplification occurring on scales up to and includ- 
ing the full star-forming disk. Assuming equipartition between the 
turbulent kinetic and magnetic energies, we estimate an upper limit of 
the equipartition turbulent B field strength of 514 uG (Methods). This is 
comparable to the estimated turbulent B field strength of 305 + 15 uG 
within the central kiloparsec of the starburst region of M82 also using 
far-infrared polarimetric observations”. This indicates that the star- 
burst activity of 9i09 could be be driving the amplification of B fields 
across the disk. 

Feedback-induced turbulence is a route to accelerating the growth 
of the seed fields, but to produce the ordered field on the kpc-scales 
observed requires a mean-field dynamo”. This mean-field dynamo 
can be achieved through the rapid differential rotation of the gas 
disk, and this provides a mechanism for the ordering of an amplified 
B field driven by star formation and stellar feedback processes. 9109 
is turbulent, intensely star-forming and rapidly rotating (maximum 
velocity Uma, 300 km s”). This suggests that rather than an episode 
of violent feedback priming a large-scale but turbulent field that 
later evolves into an ordered field during a period of relative quies- 
cence”, the small-scale and mean-field dynamo mechanisms oper- 
ate in tandem. We estimate that the mean-field dynamo in 9i09 has 
not yet had time to maintain or amplify the B field (Methods). This 
indicates that the intense starburst is most important in amplify- 
ing the galactic field at z = 2.6. We postulate that this ‘dual dynamo’ 
might be the common mode by which galactic-scale ordered mag- 
netic fields are established in young gas-rich, turbulent galaxies inthe 
early Universe. 

Coherent magnetic fields consistent with the mean-field dynamo 
have been observed at z = 0.4 by means of Faraday rotation of a back- 
ground polarized radio source? (note that such observations are not 
possible for 9i09). Magnetic fields are already known to be present 
inthe environment around normal galaxies at z = 1 as revealed by the 
association of Mg II absorption systems along quasar sightlines that 
show Faraday rotation”, and indirectly through the existence of radio 
synchrotron emission from star-forming galaxies. However, mapping 
the B fields in individual galaxies at high redshift has, so far, proved 
challenging. Our observations show that the polarized emission from 
magnetically aligned dust grains is a powerful tool to trace the B fields 
of the cold and dense ISM in high-redshift galaxies. 

Galaxy 9i09 is a particularly luminous example of a population of 
dusty star-forming galaxies in the early Universe that contribute to 
the cosmic infrared background (CIB). If the 1% level of polarization 
detected in 9io9 is representative of the general population of dusty 
star-forming galaxies” then routine detection and mapping of mag- 
netic fields in galaxies at high redshift is feasible (thatis, inintegration 
times of less than 24 h) even in unlensed systems with ALMA. This offers 
anew window to characterize the physical conditions of the ISM in 
galaxies when galaxy growth was at its maximum, and will enable a 
better understanding the role of magnetic fields in shaping the early 
stages of galaxy evolution. The strength of the galactic magnetic field 
in local spiral galaxies is of order 10 uG (ref. 1), and up to an order of 
magnitude higher in starbursts®. Without resolving the polarization 
field in 9i09 below 100 pc scales it is not possible to reliably estimate 
the B field strength using dust polarization observations. Nevertheless, 


given the injection of kinetic turbulence driven by stellar feedback we 
estimate the strength of the B field in 9i09 to be probably greater than 
that of local spiral galaxies, but similar to that of the central regions of 
nearby starburst galaxies (Methods). 

Finally, these observations imply the CIB itself may be weakly polar- 
ized”, Although misalignments of galaxies along the line of sight will 
serve to reduce the net polarization of the CIB, ifthe orientation of disks 
that host galactic-scale ordered B fields is correlated on large scales 
due to tidal alignments”, then a polarization signal could remain and 
therefore fluctuations in the polarization intensity of the CIB could be 
used as anew probe of the physics of structure formation”. This has 
consequences for cosmological experiments that seek to derive infor- 
mation on primordial conditions from observations of the polarization 
of the cosmic microwave background, especially ifa curl componentis 
present in the CIB polarization field”. A polarized component of the 
CIB at millimetre wavelengths, of extragalactic origin, dominated by 
emission atz = 2 and witha power spectrum that is shaped by large-scale 
structure at this epoch, will be a subtle but important foreground for 
future precision cosmic microwave background experiments to con- 
tend with. 
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Methods 


Observations and data reduction 

Galaxy 9i09 at celestial coordinates a= 02h 09 min 41.3 s, ô = 00° 
15’ 58.5” (J2000) was observed by the ALMA 12 m array in project 
2021.1.01461.S. The target was observed in two sessions, in which each 
session consisted of observations of sufficient duration to measure 
the rotation of the parallactic angle of the telescope by 60° or more. 
The first session, comprising two Execution Blocks, was executed on 
12 April 2022, and the second, comprising one Execution Block, was 
executed on 14 April 2022. The ALMA grid source J2253 + 1608 was 
observed as a bandpass and flux calibrator, JO006-0623 was used as the 
polarization calibrator, JO208-0047 was used as the phase calibrator 
for the first session and JO217 + 0144 was used as the phase calibrator 
for the second session. 

The data were processed inthe Common Analysis Software Applica- 
tions (CASA) v.6.2.1 using both the standard ALMA calibration pipeline 
and an extra polarization calibration script. The standard calibration 
pipeline first applied a series of steps that included corrections to the 
amplitudes (based on system temperature measurements) and phases 
(based on water vapour radiometer measurements). This was followed 
by calibration steps in which corrections for the amplitudes and phases 
as a function of frequency were derived using the bandpass calibrator, 
scaling factors for the amplitudes were derived using the flux calibra- 
tor and corrections for the phases and amplitudes versus time were 
derived using the phase calibrator. The polarization calibration script 
derived and applied further corrections for instrument-related polari- 
zation effects related to imperfections in the feeds and the orientation 
of the feeds as a function of time. The uncertainty in the final linear 
polarization calibration fraction was 0.1%, whereas the uncertainty in 
the polarization angle was 1°. 

We also used CASA for imaging of the calibrated measurement set. 
The procedure tclean was used to produce the /QUV Stokes images, 
with the ‘clarkstokes’ deconvolver, natural weighting of the visibilities 
and a common restoring beam (ensuring all Stokes images share 
the same synthetic beam). Dirty images were cleaned to a stopping 
threshold of 10 Jy and the final beam shape had a full-width at half- 
maximum of 1.2” x 0.9” at position angle 68°. All images were primary 
beam-corrected to account for fall-off in sensitivity away from the phase 
centre, but note that the target was compact and this represented a 
negligible correction to the measured flux densities. ALMA provided 
asystematic uncertainty in linear polarization of 0.03% with minimum 
detectable polarization of 0.1%, and therefore all quoted uncertainties 
in this work have the minimum detectable polarization of 0.1% added 
in quadrature. 


Measuring polarization 
To account for the vector quantity of the polarization measurements, 
we estimate the integrated polarization fraction as 


\J<Q} + (Uy? -b () 


P= 
dd) 
and the B field orientation as 
1 U Tl 
HF paretan( | + > (2) 


where <i, (Q) and (U} are the mean of the Stokes /, Q and U for measure- 
ments witha polarized intensity signal to noise ratio P/op 2 3, and bias b 
isthe standard error of the uncertainties of the Stokes Qand Uimages, 
oand oy. The mean polarization fraction, <P), and B field orientation, 
(Xs, are estimated as the mean distribution of the individual measure- 
ments of Pand x, in independent beams (that is, using the individual 
Stokes Qand Uand then computing Pand x,). The uncertainties in these 


values are estimated as the standard deviation of each distribution, 
respectively. Extended Data Fig. 1 shows the distribution of individual 
measurements at the Nyquist sampling for pixels with PI/0,, > 3. 


Lens modelling 

Our goal is to obtain the configuration of the B field morphology in 
the source plane using previous knowledge of the lens model of 9109. 
The lens model is computed using a singular isothermal ellipsoid and 
a shear component. The lens is described by the critical radius (that 
is, Einstein radius), 6,, lens offset position from the central coordinates, 
Xo Yo ellipticity, €, shear, y, and shear orientation, 6,. Extended Data 
Table 1 summarizes the best-fit lens model parameters from ref. 5. The 
morphology of the continuum thermal emission in the source plane 
is assumed to match that of the CO(4-3) molecular gas emission’. The 
CO(4-3) spectral cube can be fit with a kinematic model based on the 
ring or disk morphologies of the molecular gas reservoirs of local 
ultraluminous galaxies”, with the source plane spectral cube well 
modelled with a disk of maximum radius Rmax =322 mas 
(2, 647788, pc) and maximum rotation speed of Vnax = 300 km s”). The 
disk has a tilt angle of 5 + 4° (east of north) and is inclined along the 
line of sight by 50°23’. 

We produce synthetic polarimetric observations of our ALMA pola- 
rimetric observations as follows. First, we compute the singular iso- 
thermal ellipsoid and shear lens model from the parameters given in 
Extended Data Table 1 using LENSTRONOMY”. Using a grid of 
801 x 801 pixels? with a scale of 0.01’ pixel, which corresponds toa 
spatial scale of 82 pc in the source plane, we model the source plane 
using an asymmetric 2D Gaussian profile with a full-width at 
half-maximum equal to the Rmax and Rmin at an angle of 0 (Extended 
Data Table1). The source plane and image plane at the native resolution 
of 0.01 is shown in Fig. 2, also showing the caustics of the lens model. 
To mimic the observed data, we use the simobserve task of CASA, which 
simulates the observation of a given sky model (that is, the modelled 
image plane) with the ALMA 12 m array. We match the parameters of 
the real observations as closely as possible: we ‘observe’ in two sessions 
using antenna configuration ‘2’ for Cycle 8. The integrations start at 
hour angles of -2.339 and -2.422 h for the two sessions and last 1.86 
and 0.86 h on 12 and 14 April 2022, respectively, ensuring an identical 
sampling of the uv plane. A model is applied to simulate noise in the 
simulated observations, dominated by a thermal component, using 
the atmospheric transmission at microwaves model to simulate the 
atmospheric profile at the ALMA site, with the precipitable water 
vapour column as a scaling parameter. We assume a precipitable water 
vapour of 0.6 mm, equivalent to the average column across the obser- 
vations. We clean the simulated visibilities in the same manner as the 
real data, but apply a restoring beam equivalent to the common beam 
derived for the real observations, to ensure that the final synthesized 
beams of the simulated images match the real data. 

The B field orientation in the source plane (sp) is assumed to be con- 
stant, Xz sp with a constant polarization fraction of unity, P,,. The model 
Stokes Q,, and U,, inthe source plane are computed as 


Qp = RpCOs(2Xp sp) * lsp (3) 


Usp = Ppsin(2Xg sp) * fsp (4) 


where /,, is the Stokes /in the source plane. We multiply by /,, to convert 
the Stokes Q,, and Uin surface brightness, which allows us to compute 
the lens model. Figure 2 shows the B field orientation inthe source plane 
with vector line lengths proportional to the polarization fraction, that 
is, P,, =1%. As the Stokes Q,, and U, are density profiles, these images 
can also be lensed using the same procedure as the Stokes /,,. 

The polarized intensity, PI,,, polarization fraction, P, and B field 
orientation, Xip of the final synthetic polarimetric observation in the 
image plane are computed as 
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Pl, = Qp +U (5) 
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er (6) 
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The final synthetic observations Stokes /ip, Qip and U,, are shown in 
Fig. 1. This figure also shows the synthetic polarized intensity, PI,,,and 
the B field orientation, x,,, with the length of the lines proportional to 
the synthetic polarization fraction P,,. Note that the B field orientation 
of the synthetic observations does not have the added rotation of 90° 
as shown in equation (2). This is because we are modelling the B field 
rather than the E vector measured by the ALMA polarimetric observa- 
tions. However, to compare with our observations, Fig. 1 shows the 
synthetic Stokes Qand Uin the same reference (thatis, E vector) as the 
ALMA polarimetric observations. 


Constraining the B field orientation 

It has been shown that for a non-rotating lens the polarization vector 
of the electromagnetic wave does not rotate. This result has been dem- 
onstrated using pure geometrical definitions”, and from first princi- 
ples using a Newtonian potential and solving the spacetime metrics”. 
In general, the photons travel along null geodesics, which indicates 
that the vector properties are time invariant across the geodesic. This 
property allows us to study the intrinsic B field geometry of the source 
amplified by the gravitational lens. 

We explore potential changes of the B field orientation in our lens 
model. With a polarization fraction of 1% we set a range of constant 
Xesp Orientations in the source plane. Extended Data Fig. 2 shows the 
source plane and lensed model for several values of Xs sp = 5, 10, 45, 90° 
at a constant P,, = 1% using the same lens model as described above, 
and simulating the observations in the same way to produce convolved 
and noisy images at the same scale as the data. We conclude that the 
lens model and convolution do not produce rotation or change in 
the polarization fraction from the source plane to the image plane. 
For Xpsp = [5, 10]°, the final B, has an orientation similar to the beam 
position angle in regions with a low signal to noise ratio. However, 
the median B field orientation in the final synthetic observation is 
consistent with the modelled B field orientation in the source plane, 
with the presence of noise and asymmetric synthetic beam contrib- 
uting to the angular dispersion in the image plane. We find that the 
background noise and uv plane sampling (in general, imaging of the 
visibilities) introduce an uncertainty of roughly 10° in <xz;,), in which 
the uncertainty is dominated by the signal to noise ratio of the ALMA 
polarimetric observations. 

The B field orientation can also be constrained using the Stokes Q 
and Uimages. Our observations show that Stokes Q is negative and 
Stokes Uis consistent with zero (Fig. 1). This result shows that the E 
vector is mainly inthe north-south direction. This configuration gives 
us the opportunity to tightly constrain the possible range of values in 
the B field orientation. Extended Data Fig. 3 shows the Stokes /, Q and 
Uimages, and the polarized intensity of synthetic observations for the 
same set of artificial B field orientations in the source plane. Note the 
change from negative to positive Q from x,,, of 5 to 90°. This behaviour 
is expected, but we show it for completeness. We can constrain the 
orientation by rotating the B field until the corresponding Stokes U 
image shows a 3o detection in the synthetic polarimetric observations. 
This condition is met when x,,,, = + 10°. For orientations deviating out- 
side +10° the synthetic observations are not consistent with our obser- 
vations. We conclude that the most likely B field orientation consistent 


with our observations is X, ,, = 5"j, where we have set the B field ori- 
entation to be parallel to the CO(4-3) emission of the disk inthe source 
plane’. 


Constraining the polarization fraction 

We assume a constant polarization fraction of 1% in the source plane. 
The lens model shows that the median polarization fraction is con- 
sistent with this level without dispersion (Extended Data Fig. 2). We 
conclude that the lens model does not change the polarization frac- 
tion from the source plane to the image plane. Using the simulated 
observations, we estimate the median and r.m.s. polarization fraction, 
finding that the combination of uv-plane sampling and background 
noise adds an uncertainty of roughly 0.1-0.2% in the polarization 
fraction. 


Energy equipartition 

We estimate the B field strength assuming equipartition between 
the turbulent kinetic energy and the turbulent magnetic energy. Let 
the turbulent kinetic energy, U,, and turbulent magnetic energy, 
Us, be 


1 

U= 59% (8) 
2 

Uy= a (9) 


where p is the volume density, o, is the velocity dispersion and Bis the 
field strength. Then, assuming equipartition, U, = Ub, the field strength is 
Beg = [4119 0. (10) 

To estimate the baryonic volume density, p = M/V, we use the disk 
volume, V= 1h, and mass of the molecular gas of the galaxy. The 
radius isr=2, 647788, pc and the height is h < 600 pc (ref. 5); note that 
the disk height is unresolved in the observations of the CO molecular 
gas, so we take the size of the beam of the observations as a conserva- 
tive upper limit, also noting that the scale heights of bursty disk galax- 
ies at high redshift are observed to be larger than local disks, with 
heights of several hundred parsecs typical. The total molecular mass 
is estimated from the CO luminosity’ as My, = (7.540.1) x 10!° M, and 
the velocity dispersion of the gas is 0, = 73 + 4 kms‘, from kinematic 
fitting of the disk model to the CO data cube’. From this, we estimate 
an equipartition B field strength of B,, < 514 uG, however note that the 
observed velocity dispersion may be affected by large-scale flows from 
galactic winds and shearing effects by the rotating galactic disk, such 
that the turbulent component velocity dispersion might be lower. 
These effects produce an overestimated measurement of the B field 
strength”. 

Spiral galaxies have an average ordered B field strength of around 
5+2 uG with a total B field strength of 17 + 14 uG assuming equiparti- 
tion between the total B field and total cosmic-ray electron density’. 
Arevised equipartition formula to account for energy losses in nearby 
(within 160 Mpc) starburst galaxies estimated” equipartition B field 
strengths in the range 70-770 uG. Recent far-infrared polarimetric 
observations” of the starburst galaxy M82 showed that the turbulent 
Bfieldstrengthis 305 + 15 uG within the central kiloparsec region. These 
authors modified the Davis-Chandrasekhar-Fermi method to account 
for the large-scale flow of the galactic outflow. The correction from the 
equipartition B field strength, Beq = 540 + 170 uG, was estimated to be 
roughly 25%. 


Ordering timescale 
We estimate the timescale to order a large-scale magnetic field ina 
disk. The ordering timescale can be estimated as” 


Pee ee 

te Bh Pd 
where his the half-thickness of the galactic disk and l, is the large-scale 
coherence length. £ is the turbulent diffusivity defined as £ = lo,/3, 
where lis the coherence length of the small-scale turbulence and o, is 
the turbulent velocity. D, is the dynamo number defined as 
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where Qis the angular velocity. 

For 9i09 we assume a Gaussian scale height? h = 300 pc and 
0,=73+4kms ‘ofthe molecular gas. To estimate the angular velocity 
we use the deprojected circular velocity, Uma,=360' kms", at 
the maximum radius of the disk, fnax=2, 64778 pc, yielding 
Q = Uimax/ Tmax = 145+8kms kpc”. 

The typical coherence length of the small-scale dynamo is driven 
by stellar activity in the molecular gas of galaxies on scales of 
l=1-10 pc (ref. 36). We assume a large-scale coherence length in the 
range l, = 0.5 - 2.0 kpc, corresponding to 50-2,000 times larger than 
the turbulence coherence length, /. We estimate a dynamo number 
D4=3.2 + 0.4 and turbulent diffusivity £ = (4.2 + 2.0) x 10” cm? s”. 
Finally, the ordering timescale is estimated to be in the range of 
tx 0.4-18 Gyr, noting that only 2.5 Gyr have elapsed by z = 2.6. Analyti- 
cal solutions of the evolution of B fields in spiral galaxies predicts that 
a large-scale B field can be formed within 2-5 Gyr (ref. 37). This study 
suggests that the mean-field dynamo may be already active at z < 3. 
Note that 9i09 has a very high star-formation rate, which may inject 
higher turbulent energy into the system than those studied in ref. 37 
affecting the ordering of the B field in the galaxy’s disk. The regular B 
field is amplified if the dynamo number is larger than the critical value 
of Dy, 7 estimated from numerical simulations of galactic dynamo 
models”. For comparison, the Milky Way has a dynamo number of 
Dg=9 > Dac Which shows that large-scale dynamo mechanism is 
important in our galaxy. For 9109 we estimate D, = 3.2 + 0.4, below 
the critical value. In addition, the time for 9i09 to complete a rotation 
is t= 45 Myr, which indicates that roughly 9-400 galactic rotations 
at a galactocentric radius of 2.6 kpc are required to order the large- 
scale B field. 


Data availability 


The ALMA data used in this work will be available through the ALMA 
Science Archive (https://almascience.nrao.edu/aq) with reference to 
the project code 2021.1.01461.S. 
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Extended Data Fig. 1| Histograms of the measured polarization fraction individual measurements are taken in independent beams. The median and lo 
and magnetic field orientation. The polarization fraction (a) is displayed in scatter of the distributions is given for each measurement. 


bins of 0.5% and the B field orientation (b) in bins of 10 degrees, where 
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Extended Data Fig. 2 | Characterisation of the magnetic field orientation caustics in the source plane (green line), and image plane (yellow line). The 
and polarization fraction. The B field orientation (white lines) at 5 degrees (a), median andr.m.s. of the polarization fraction and B field orientation are shown 
10 degrees (b), 45 degrees (c), and 90 degrees (d) are shown over the source at the top right of each panel. The red ellipse indicates the synthetic beam of 


plane (a), lens model (e), and synthetic observations (f). We indicate the the simulated observations. 
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Extended Data Fig. 3 | Variation of the Stokes Q and U parameters with 
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polarimetric observations. The red ellipse shows the synthetic beam of the 


simulated observations (a). For Stokes /, the contours increase as g, x 2°45, For 
Stokes Qand U, and PI, the contours start at 3g and increase in steps of lo. Note 
that at o,, = 10 degrees a 30 signal becomes detectable in Stokes U. 


magnetic field orientation. The B field orientation (white lines) at 5 degrees 
(a), 10 degrees (b), 45 degrees (c), 90 degrees (d) are shown over the Stokes / (a), 
Stokes Q (e), Stokes U (f), and polarized intensity (g) of the synthetic 


Extended Data Table 1| Lens model parameters 


Einstein radius Or 2.45” 
Ellipticity € 0.072 
Ellipticity angle Oe 67.2° 
Shear y 0.056 
Shear orientation 0, —91° 
Right Ascension x, 02°09"41.2545 
Declination Ye +00°1558.450” 


Angles are measured counterclockwise, i.e., East of North. Coordinates in J2000. 


